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5 Efficient tests of association for quantitative traits and affected-unaffected 

studies using pooled dna 

Related Application 

This application claims priority to U. S. Ser. No. 60/238,381, filed October 6, 2000 

[21402-139] which is incorporated herein by reference in its entirety. 

1 0 Background of the Invention 

The complex diseases that present the greatest challenge to modern medicine, 
including cancer, cardiovascular disease, and metabolic disorders, arise through the interplay 
of numerous genetic and environmental factors. One of the primary goals of the human 
genome project is to assist in the risk-assessment, prevention, detection, and treatment of these 

1 5 complex disorders by identifying the genetic components. Disentangling the genetic and 
environmental factors requires carefully designed studies. One approach is to study highly 
homogenous populations (Nillson and Rose 1999; Rabinow, 1999; Frank 2000). A recognized 
drawback of this approach, however, is that disease-associated markers or causative alleles 
found in an isolated population might not be relevant for a larger population. An attractive 

20 alternative is to use well-matched affected-unaffected studies of a more diverse population 

Even with a well-matched sample set, the genetic factors contributing to an aberrant 
phenotype may be difficult to determine. Traditional linkage analysis methods identify 
physical regions of DNA whose inheritance pattern correlates with the inheritance of a 

25 particular trait (Liu 1997; Sham 1997, Ott 1999). These regions may contain millions of 

nucleotides and tens to hundreds of genes, and identifying the causative mutation or a tightly 
linked marker is still a challenge. A more recent approach is to use a sufficiently dense 
marker set to identify causative changes directly. Single nucleotide polymorphisms, or SNPs, 
can provide such a marker set (Cargill et al. 1999). These are typically bi-allelic markers with 

30 linkage disequilibrium extending an estimated 10,000 to 100,000 nucleotides in heterogeneous 
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can provide such a marker set (Cargill et al. 1999). These are typically bi-allelic markers with 
linkage disequilibrium extending an estimated 10,000 to 100,000 nucleotides in heterogeneous 
human populations (Kruglyak 1999; Collins et al. 2000). Tens to hundreds of thousands of 
these closely spaced markers are required for a complete scan of the 3 billion nucleotides in 
5 the human genome. Because each SNP constitutes a separate test, the significance threshold 
must be adjusted for multiple hypotheses (p-value ~ 10" 8 ) to identify statistically meaningful 
associations. Consequently, hundreds to thousands of individuals are required for association 
studies (Risch and Merikangas 1996). 

10 The most powerful tests of association require that each individual be genotyped for 

every marker (Fulker et al. 1995, Kruglyak and Lander 1995, Abecasis et al. 2000, Cardon 
2000) and remain far too costly for all but testing candidate genes. An alternative that 
circumvents the need for individual genotypes, related to previous DNA pooling methods for 
determination of linkage between a molecular marker and a quantitative trait locus (Darvasi 

15 and Soller 1994), is to determine allele frequencies for sub-populations pooled on the basis of 
a qualitative phenotype. Populations of unrelated individuals, separated into affected and 
unaffected pools, have greater power than related populations. Limited guidance has been 
provided, however, regarding the sample size requirement of tests using pooled DNA relative 
to individual genotyping, or the efficiency of tests based on a quantitative phenotype relative 

20 to an affected/unaffected design. 

The phenotypes relevant for complex disease are often quantitative, however, and 
converting a quantitative score to a qualitative classification represents a loss of information 
that can reduce the power of an association study. The location of the dividing line for 
25 affected versus unaffected classification, for example, can affect the power to detect 

association. Furthermore, pooling designs based on a comparison of numerical scores are not 
even possible with a qualitative classification scheme. These distinctions can be especially 
relevant when populations contain related individuals and qualitative tests have a disadvantage 
(Risch and Teng 1998). 

30 

When performing risk assessment to determine whether a person suffers from or is at 

risk of developing a complex disorder often requires measuring an underlying quantitative 

phenotype. Association studies in unrelated populations can implicate genetic factors 
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contributing to disease risk, and experiments using pooled DNA provide a less costly but 
necessarily less powerful alternative to methods based on individual genotyping. Association 
studies require markers in linkage disequilibrium with causative genetic polymorphisms. 
Although the sample sizes required for pooling and individual genotyping studies have been 
5 compared in certain instances, general results have not been reported in the context of 

association studies, nor have there been clear comparisons of pooling based on quantitative 
and qualitative (affected/unaffected) phenotypes. Association tests of DNA pooled on the 
basis of a quantitative phenotype are analogous to selection experiments for quantitative trait 
locus (QTL) mapping. For a QTL with a weak effect on a phenotype, the mean phenotypic 
10 value of individuals selected to exceed a threshold is proportional to the mean allele 

enrichment. This suggests that genotyping of a certain percentage of the upper and lower 
phenotypic values of an unrelated population is useful to estimate the effect of a marker on a 
quantitative phenotype, such as in pooling studies. There is a need in the art to examine the 
sample size requirements of association tests for quantitative traits using pooled DNA. 

1 5 Summary of the Invention 

The present invention is based, in part, on the discovery of methods to detect an 
association in a population of individuals between a genetic locus and a quantitative 
phenotype, where two or more alleles occur at a given genetic locus, and the phenotype is 
expressed using a numerical phenotypic value whose range falls within a first numerical limit 

20 and a second numerical limit. These limits are used to provide for subpopulations that consist 
of upper and lower pools. 

In some embodiments, the population of individuals includes individuals who may be 
classified into classes. In certain aspects of the invention, these classes are based on age, 
gender, race, or ethnic origin. In other aspects, some or all members of a class are included in 

25 the pools. 

In various embodiments, these numerical limits are chosen so that the upper pool 
includes the highest 19%, 27%, or 37% of the population. In other embodiments, the 
numerical limits are chosen such that the lower pool includes the lowest 19%, 27%, or 37% of 
the population. 

30 In some embodiments, the upper and lower pools have the same number of individuals. 
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In one embodiment of the invention, the numerical limits are chosen to correlate with 
error of measurement determinations. In some embodiments, the numerical limit on the error 
of measurement is about 0.04 or about 0.01. 

In some embodiments, methods to detect an association in a population of individuals 
5 between a genetic locus and a quantitative phenotype are useful to determine the genetic basis 
of disease predisposition. 

In other embodiments, the genetic locus analyzed contains a single nucleotide 
polymorphism. 

In the present invention, the population of individuals can include unrelated 
10 individuals. 

Unless otherwise defined, all technical and scientific terms used herein have the same 
meaning as commonly understood by one of ordinary skill in the art to which this invention 
belongs. Although methods and materials similar or equivalent to those described herein can 
be used in the practice or testing of the present invention, suitable methods and materials are 
1 5 described below. All publications, patent applications, patents, and other references 

mentioned herein are incorporated by reference in their entirety. In the case of conflict, the 
present specification, including definitions, will control. In addition, the materials, methods, 
and examples are illustrative only and not intended to be limiting. 

Other features and advantages of the invention will be apparent from the following 

20 detailed description and claims. 

Brief Description of the Figures 

Fig. 1 . The sample size required to achieve a type I error rate of 5xl0" 8 and a power of 0.8 for 
a QTL for a complex trait is shown for pooled DNA designs relative to individual genotyping. 
The ratio N c< /N' m &\ y for affected-unaffected pools (dashed line) is shown as a function the 
25 disease incidence r, while the ratio Akij/Windiv (solid line) is shown as a function of the fraction 
p of the total population selected for each pool. The optimum value of Maii/7Vi nc iiv is 1 .24, 
occurring at p= 27% selected for each pool. 
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Fig. 2a Exact numerical results for the sample size N required to achieve a type I error rate of 
5xl0" 8 with a power of 0.8 are shown for affected-unaffected pools (dashed line) and tail pools 
(solid line) as a function of the additive variance, or equivalently the genotype relative risk for 
5 a heterozygote, for an allele with frequency 0.1 and purely additive variance. Analytic 

approximations (solid circles), Eqs. 1 and 2, are indistinguishable from the exact results when 
the genotype relative risk is smaller than a factor of 2. The disease incidence r is 10% for the 
affected-unaffected pools, and 27% of the population is selected for the each of the tail pools. 

10 Fig. 2b The frequency difference at the significance threshold is shown for the same 

parameters as panel a. This threshold determines the measurement accuracy required for an 
association test based on pooled DNA. 

Detailed Description of the Invention 

1 5 The present invention provides analytic results for association tests. It is shown that the 

results obtained closely approximate the analytic results to exact numerical calculations. The 
invention further extends the analysis to qualitative phenotypes using a genotype relative risk 
model. 

20 A particular quantitative phenotype Xis standardized to have unit variance and zero mean. 

The phenotype is hypothesized to be affected by alleles A\ and A 2 , with frequencies p and \-p 

respectively, at a particular QTL. The population fractions P(G) for genotypes G = A\A\, 

A\Ai, and A2A2 are assumed obey Hardy-Weinberg equilibrium. Using standard notation for a 

variance components model (Falconer and MacKay, 1996), the effect juc of genotype G on 
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phenotype X is a-p for A \A \ , for A \A 29 and -a-// for A2A2. The constant // = (2p-\ )a 
+ 2p(l-p)d ensures that the mean of X is zero. The ratio dla describes the inheritance mode 
for allele A\. Dominant, recessive, and additive inheritance are special cases with dla equal to 
+1,-1, and 0, respectively. 

5 

The phenotypic variance due to the QTL may be partitioned into the additive variance <r A 2 and 

the dominance variance ab 2 , with 

cta 2 + <7 D 2 = 2pq[a^d(p^q)] 2 + 4p 2 q 2 d l . 

The additive variance is often much larger than the dominance variance even if the inheritance 
1 0 mode is not purely additive. The exceptions are QTLs with a recessive minor alleles and 
dominant major alleles, which are difficult to detect in unselected populations. The 
contribution of remaining genetic and environmental factors is assumed to follow a normal 
distribution with residual variance or 2 , 
or 2 = H<7 A 2 +ob 2 ). 

1 5 Of particular interest here are complex traits: the effect of any single QTL is small, cr A 2 +ot> 2 < 
0.05, and the residual variance <r R 2 is nearly 1. 

A genotype relative risk model corresponds to classifying individuals as affected {X> Xj) or 
unaffected (X< X T ) based on a specific threshold Jf/-. The proportion r of the total population 
20 that is affected is the overall risk or disease incidence; the probability that an individual with 
genotype G is affected, relative to the probability for an individual with genotype A 2 A 2i is the 
genotype relative risk. If the inheritance mode of A\ is additive and a is small compared to or, 
the relative risk is multiplicative with allele dose. 
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The sample size N required to detect association between genotype G and the quantitative 
phenotype or the disease risk depends on the type I error rate a, the type II error rate /?, and the 
test statistic and experimental design (Snedecor and Cochran, 1989), as well as on the 
underlying genetic model For a one-sided test of a single marker, a = 1 - <P(z a ), where <D(z) 
5 is the cumulative probability distribution for standard normal deviate z, defines a in terms of 
deviate z a . Similarly, l-/?is the power to reject the null hypothesis and z\-p = <£ _1 (/?). For a 
genome scan, the values a= 5xl(T 8 (z a = 5.33) and 1-/?= 0.8 (zi_/? = -0.84) have been 
suggested (Risch and Merikangas, 1996). 

10 We consider two experimental designs using DNA pooled from individuals selected from a 
sample of size N: affected-unaffected pools, with DNA pooled from n affected and n 
unaffected individuals; and tail pools, with DNA pooled from n individuals at each tail of the 
phenotype distribution. The test statistic for these designs is the frequency difference of the A\ 
allele between the pools. The multinomial distribution describing the test statistic may be used 

1 5 to calculate exactly the sample size required to achieve statistical significance at specified 
power. 

When the number of A\ alleles summed over both pools is large, the distribution of the test 
statistic is approximately normal. A significant association is detected if the allele frequency 
20 difference between pools is at least z a times the standard deviation of its estimator, or z a p 1/2 (l- 
p) m ln m . Furthermore, when the additive variance a A 2 is small and the residual variance 
or 2 is close to 1, convenient analytic approximations for the sample size requirements may be 
derived. 
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For the affected-unaffected design, n = rN of the individuals are expected to be diagnosed as 
affected, and an additional n matched controls are selected from the remainder of the 
population. The analytic approximation for the sample size is 

N„ = [z a -z x . p f [okW] • 2/<l-r)V[l +^ 7 <l-a- R 2 ) ,/2 /2 3/2 o- R 2 / > 1/2 (l-p) 1/2 ] 2 . (Eq. 1) 
5 The termy is the height of the standard normal distribution at the normal deviate XjJok 
corresponding to the threshold between affected and unaffected phenotypic values. 

The tail pools are parameterized by the fraction p = nIN of population selected for each pool, 
and p plays a role analogous to the overall disease incidence r in the affected-unaffected 

10 design. The analytical approximation for the sample size is 

Mail = [z a - zx-pf [<7r 2 /<Ta 2 ] • p/2y\ (Eq. 2) 

where y is the height of the standard normal distribution for normal deviate <&~ ] (p). The 
design may be optimized by selecting p to minimize Mail, which corresponds to minimizing 
p/2y . With this approximation, the optimal fraction is 0.27 and is independent of a, /?, and all 

1 5 parameters of the genetic model. 

A third method, individual genotyping, serves as a baseline for evaluating the efficiency of the 
two pooling-based methods. The sample size required to achieve significance using individual 
genotyping is 

20 JVindiv = [z a - zi^r] 2 /ox 2 , (Eq. 3) 

based on a regression model of phenotypic value on allele dose. 



Detailed Description of Analytical Methods 

The genotype-dependent phenotype distribution in the variance components model is 
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P(X\G) = (2^r I/2 expK^) 2 /aR 2 ], 

and the overall phenotype distribution is the sum of the three normal distributions, 
P(X) = X G P(X\G)P(G). 

When an upper threshold X\j is specified to select a fraction p of the total population with 
5 phenotypic values above the threshold, the equation 
p = Z G {1 -<S>[(Xu-MG)/crG]}P(G). 

may be solved numerically for X v as a function of p. The genotypes of individuals selected by 
X>Xy follow a multinomial distribution; the probability that an individual has genotype G is 
6u(G) ={ \-Q>[(Xu~MG)lcrG]}P{G)lp. A multinomial distribution is similarly defined using a 
1 0 lower threshold X L , 

1 = E G 9l{G) = p'lc ®[(X l -iig)/<t g ]P(G). 

For an affected-unaffected design, the fraction in the upper pool is r and the fraction in the 
lower pool is 1-r, yielding X v = X L = X T . The relative risk for genotype G is [OiAG)IP(G)\l 
[e l AA 1 A 2 )IP(A 1 A 1 )\ 

15 

Sample size requirements may be obtained directly from the multinomial distributions of 
genotypes by exhaustively tabulating allele counts C v and Q in the upper and lower pools for 
each distinct composition of genotypes among the n selected individuals. The distribution 
corresponding to null hypothesis, 6(G) = P(G), is used to define the smallest threshold AC 
20 such that Cu~C L > AC with probability a or less. The discrete allele count usually yields the 
strict inequality. Next, the distributions under the alternative hypothesis are considered, and 
the probability that C\j-Cl> AC is tabulated to provide the power. If the power is greater 
than or equal to the specified 1-J3, the choice of n and Af = nip or nlr is feasible. A search is 
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performed for the smallest feasible N with r or p specified. For tail pools, p is then varied to 
find the overall optimum. 

When the number of alleles summed over both pools is large, the allele frequency difference 
5 follows a normal distribution. Under the null hypothesis, the mean is zero and variance is 
oo 2 /n =p(\-p)/n. This result is derived by noting that the variance of the frequency difference 
is twice the variance of the mean for a single pool of n individuals. The allele frequency 
variance for an individual is p(\~p)/2, and averaging over the n individuals reduces the 
variance by the factor «. Under the alternative hypothesis, the expected allele frequency 
10 difference zip is 

4> =Pu-Pl = ^g[0iAG)~0l{G)]pg 

where the genotype-dependent allele frequency pc is 1 for G = A\A\ 9 0.5 forA\A2, and 0 for 
A2A2. The variance is <j\ 2 /n 9 where o\ is obtained from the multinomial distribution (Beyer, 
1984), 

15 <r, 2 = S c [0a(G) + 6 l (G)]pg - (pa 2 +p L 2 ). 

The number of individuals required per pool for type I error a and power 1-/3 is 
n = [z a a 0 -zi-p<Ji] 2 /Ap 2 . 

For affected-unaffected pools, N = nlr is the required sample size. For tail pools, N = nip, and 
p is varied to find the smallest N. 

20 

The normal approximation underestimates the sample size requirement relative to the exact 
results from the multinomial distribution. When the sum of the alleles in both pools is at least 
60, the difference in sample sizes is no greater than 5%. We chose 60 alleles in both pools as 
the criterion for switching from the multinomial to the normal calculation. Standard algorithms 
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were employed to perform the root search for X\j and X L , the optimization, and the integration 
over the tail of a normal distribution (Press, 1997). 

The analytic results are obtained by setting <j\ to <7o 2 and expanding Ap to second order in the 
5 effect size jug, corresponding loosely to a perturbation theory for probability distributions 
(Chandler, 1987). From a Taylor series expansion, 
0(z-6) = 0>(z) - by - ( 1/2)6 V, 

where y = (2;r)~ 1/2 exp(-z 2 /2). Substituting this result into the expressions for 6(G) using b = 
/^/ctr and z = XJa^ = Q>~\l-p) 9 where X is the threshold used to select the pool, yields for the 
10 tail design 

Pu=P + iylp) E[(jjg/ct k )p g ] + (y\z\/2p) E[(mg/<Tr?pg] and 
PL =P ~ iylp) E[(jig/ct r )pg] + (y\z]/2p) E[(pc/a R ) 2 p G l 

The corresponding results for the affected-unaffected pools, with z = 0 _1 (l-r), are 
Pu=P + iylr) E[(tiG/a R )p G ] + (y|z|/2r) E[(//^a R ) 2 /? G ] and 
15 />i =p - W-r)]E[{MGl(JK)PG\ - [y|z|/2(l-r)] E[(/^ R ) 2 p G ]. 
The required expectation values are 
V[McPg] = ^ G P{G)pcPG = 0a[p(1- j p)/2] ,/2 , and 

E[A; 2 = XgP(G)pg 2 Pg = (1/2)(1-ct r 2 ) - 4/> 2 (l-/>) 2 aJ + (2p-l)a D 2 /2 « <7 A 2 /2. 
The results for 4p, 
20 49 = 2 m yaoa A /p<jR , tail pools, and 

4? = [l-h\V<T A /2 3/2 c7ocr R 2 ] ;/a 0 cr A /2 1/2 r(l-r)cr R , affected-unaffected pools, 
lead directly to Eqs. 1 and 2. 
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Approximate genotype relative risks may also be obtained from the Taylor series expansion 
for 6(G). To lowest order, the relative risk for the heterozygote is approximately 1 + 
(d+a)y/roR 9 and for the A\A\ homozygote is 1 + 2ay/ra K . For additive inheritance, d = 0, and 
the relative risk is multiplicative with allele dose when ay/ra R is small. For a complex trait or 
5 is close to 1 , and for a minor allele, a « <j A /(2p) m . When the disease incidence is 10%, the 
parameter required to be small is 1 .24<7 A /p I/2 . 

For individual genotyping, the regression model used to test significance is 
X=b\(p G -p) + s, 

10 where the residual contribution € to the phenotype has zero mean and is uncorrected with pa- 
using standard statistical methods (Snedecor, 1989), the test statistic b\ under the null 
hypothesis has mean zero and variance Var(2>i|null) given by 
Var(6,|null) = AT 1 Vm{X)Nm(p g ) = \IN\p(\-p)l2\. 
Under the alternative hypothesis, the expectation for the test statistic is 

15 E(Z>0 = Cov(X, Pg )Nm{X) = a A [p(l-p)/2f\ 
and its variance is 

Var(6,|alt) = AT 1 Var(*)/Var(p G ) = <r R 2 / N b(l-/>)/2]. 

The sample size required for a one-sided test of b\ with Type I error a and power l-/?is 
N= [z a Var(6i|null) 1/2 -z^^Var^^al^^^/E^O 2 , 
20 which is the result provided in Eq. 3. 

Application of the methods of the invention 

The sample sizes required for the pooled DNA designs are compared in Fig. 1 to the sample 

size Nindiv required by individual genotyping. The ratio N c JN m & v (dashed line) is a function of 
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the disease incidence r, while MaiyWindiv (solid line) is a function of the pooling fraction p. For 
typical disease incidence, r ~ 10%, the affected-unaffected design requires a sample 5.3x 
larger than that required for individual genotyping. Compared to the tail design, it measures 
an allele frequency difference that is half as large and is approximately 4x less efficient. The 
5 tail design, with p = 27%, requires a sample only 1 .24 x larger than required for individual 
genotyping. The tail design is also robust to variation in p near its optimum, as values from 
19% to 37% drop the efficiency no more than 5%. 

The analytic theory indicates that the additive variance ox 2 , or equivalently the genotype 
10 relative risk for an allele of known frequency, is the most important factor determining the 
sample size requirements. This dependence is shown in Fig. 2a with exact numerical results 
for affected-unaffected pools (dashed line) and tail pools (solid line) for type I error of 5x1 0~ 8 
and power of 0.8. The minor allele frequency is 10%, its effect on the quantitative phenotype 
is purely additive, and the disease incidence is 10%. The analytic approximations (solid 
1 5 circles) from Eq. 1 and 2 are nearly indistinguishable from the exact results when the genotype 
relative risk drops below a factor of 2. As predicted by the analytic theory, the tail pools 
require smaller sample sizes than the affected-unaffected pools, and the gap grows wider for 
alleles with a smaller effect on the phenotype. For relative risks of 2 to 5, the deviations from 
analytic theory are moderate; above a relative risk of 5, the phenotype is monogenic with 
20 respect to locus G, and the analytic approximations for complex traits are no longer valid. 

The allele frequency difference between pools at the significance threshold is shown in Fig. 2b 
for affected-unaffected pools (dashed line) and tail pools (solid line). The measurement error 
in the allele frequency difference must be smaller than the significance threshold to detect 
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association (Darvasi, 1994). Evaluations that provide a frequency difference measurement 
accurate to 0.04 can detect association with alleles responsible for 1% of the total phenotypic 
variance, corresponding to a heterozygote relative risk of 1 .5. The allele frequency difference 
measurement must be accurate to 0.01 to detect association with an allele explaining 0.1% of 
5 the phenotypic variance, corresponding to a relative risk of 1 . 14. 

To test the range of validity of the analytic estimates for pooling, we performed a series of 
exact calculations of sample size requirements as a function of p and d/a. Large deviations 
were seen only when the magnitude of a gene effect //g approached or in size, or, 
1 0 equivalently , when <t a 2 was larger than the minor allele frequency or when a genotype relative 
risk was larger than 5 (results not shown). For additive contributions from a minor allele, the 
range of validity corresponds to cr A 2 < 2p. 



The advantages of the methods disclosed herein include the following. The optimal fraction 
1 5 for tail pooling, 27%, is independent of all model parameters including allele frequency, 

inheritance mode, effect size, and type I error and power, for virtually any QTL contributing to 
a complex trait. The exceptions to this finding are rare QTLs with relative risks of 5 or 
greater, and rare, recessive alleles, both of which are more difficult to detect than more 
frequent alleles contributing to the same overall phenotypic variance. In addition, the tail 
20 design is approximately 4-fold more efficient than the affected-unaffected design and requires 
a sample size only 24% larger than for individual genotyping. Still further, DNA pooling 
studies designed according to the present procedures disclosed herein provide extremely 
efficient methods for large-scale screening and should help to make feasible genome-wide 
association studies. 
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OTHER EMBODIMENTS 

While the invention has been described in conjunction with the detailed description 
thereof, the foregoing description is intended to illustrate and not limit the scope of the 
invention, which is defined by the scope of the appended claims. Other aspects, advantages, 
and modifications are within the scope of the following claims. 
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